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ABSTRACT 

Recent results comparing interacting galaxies to the mass-metallicity relation show that their nuclear 
oxygen abundances are unexpectedly low. We present analysis of N-body/SPH numerical simulations 
of equal-mass mergers that confirm the hypothesis that these underabundances are accounted for by 
radial inflow of low-metallicity gas from the outskirts of the two merging galaxies. The underabun- 
dances arise between first and second pericenter, and the simulated abundance dilution is in good 
agreement with observations. The simulations further predict that the radial metallicity gradients of 
the disk galaxies flatten shortly after first passage, due to radial mixing of gas. These predictions will 
be tested by future observations of the radial metallicity distributions in interacting galaxies. 
Subject headings: galaxies: abundances — galaxies: evolution — galaxies: interactions — galaxies: 
ISM 



1. INTRODUCTION 

Galaxy interactions are major actors in the process of 
galaxy evolution. Because of this leading role, the impact 
of interactions on the chemical evolution of galaxies has 
become a subject of intense scrutiny. We know little 
about heavy elements in interacting galaxies, compared 
to our understanding of more isolated systems. 

One place to look for signatures of chem ical evolution is 
the mass-metallicit y rela tion of galaxies (jLeaueux et al.l 
[19791: IRubin et alJliOSl . Low -mass galaxies contain 
fewer heavy elements than high-mass galaxies, due 
to differing star formation histories c ombined with 
the action of gas in flows and outflows (G arnettl l2002t 
iTremonti et al.l 120041 ) . In particular, star formation en- 
riches the interstellar medium via nucleosynthesis; gas in- 
flows of relatively metal-poor gas, either from outside the 
galaxy or interior to the galaxy, can change the galaxy's 
metal distribution; while outflows prevent furt her enrich- 
ment by blowing out recently produced metals ( Edmunds 
19901: lEdmunds fc Greenhowlll995l; iKoppen fc Edmunds! 
1999t lDaTcantonll2007t) . ~ 



Recent studies reveal that galaxies involved in ma- 
jor interactions fall below this fundamental relation. 
In other words, the centers of interacting galaxies are 
underabundan t com p ared to galaxies of similar mass 
(iKewlev et al.1 [20061: iRupke et all l200l lEllison et all 



120081: iMichel-Dansac et all 120081: IPeeples et al.l 120091) . 

The physical model for these deviations from the M — 
Z relation is attractively straightforward. If one or 
both of the progenitor galaxies is a spiral galaxy, ra- 
dial abundance gradi e nts exist in the mer ging galaxies 
(|Zaritskv et al.l 119941: Ivan Zee et ail 119981) . with lower 
abundances at larger galactocentric radius. During 
the merger, low-metallicity gas from the galaxy out- 
skirts is torqued into the high-metallicity galaxy ce nter 
(jMihos fc Hernquist 1996; B arnes fc Herna uist 1996|), re- 
sulting in gas with a lower average abundance. 

This gas redistribution should also result in a change 
in the gas metallicity profile across the entire disk. 
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Thus, radial metallicity gradients in mergers will differ 
from those in normal spirals. Barred galaxies, which 
are undergoing radial redistribution of gas due to bar- 
induced gas motions, indeed show evidence of cither 
shall ower gradients than are typical of unbarred spi- 
rals (iVila-Costas fc Edmundslll992l: iZaritskv et al.lll994t 
IDutil fc R.ovf Tl9991 or flatten ed oxygen gradients outside 
the c orotation radius ([Martin fc Rovl[i995l : IRov fc Walshl 
1997). They also shower lower central metallicity 
than g alaxies of similar m orphological type and lumi- 
nosity (IDutil fc Rovl 119991) . as is observed in mergers. 
iFriedli et all ([19941 ) model this behavior and demon- 
strate that the metal gradient changes are due to gas 
redistribution by the bar potential. A brief mention of 
metal redistribution in paired galaxies in the context of 
a coarse-gra i ned c osmological simulation is also made in 
iPerez etall ([20061) . However, the idea of radial gas pro- 
file changes has not been quantitatively explored using 
detailed numerical simulations in the context of strongly 
interacting galaxies, which are undergoing even more vi- 
olent gas redistributions than barred spirals. 

We here analyze N-body / smoothed-particle hydrody- 
namics (SPH) simulations of galaxy mergers to quantita- 
tively study metal redistribution due to merger-induced 
gas motions. We wish to explore nuclear metallicity and 
radial metal gradients as a function of time in the con- 
text of a merger of equal-mass spirals. This particu- 
lar scenario is relevant to the formation of ultralumi- 
nous infrared galaxies (ULIRGs) a nd QSOs, which are 
thought to arise from such a merger ([Sanders et al.lll988l : 
iVeilleux et aDl2002f ). g2] introduces our simulations and 
analysis. In § §3] — S] we present and discuss our results, 
and we summarize in $5] 

2. SIMULATIONS AND ANALYSIS 

Our starting point is eight N-body/SPH simulations 
of close-passage, equal-mass disk galaxy mergers. These 
simulations use the setup described in in §3.1 and Ap- 
pendix A of lBarnesI (|2004l ). except that ours exclude star 
formation. In brief, 87040 particles (24576 gas) were 
used. Gas mass fractions (as a fraction of the stellar+gas 
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Fig. 1. — Change in nuclear (R/R^ < 0.5) metallicity as a function of time. The metallicity is expressed in terms of Z, the mass ratio 
of an element to that of hydrogen. Strong evolution is observed between first and second pericenter, of typical magnitude 0.2 dex. The 
change in nuclear abundance is not strongly dependent on geometry, and correlates very well with the amount of gas inflow (Figure [2} • 
Vertical dashed lines locate the time of each pericenter, and data beyond the second pericenter is not plotted because of intermingling of 
the two gas disks. 



mass) were set at 12.5%, and the gas was distributed in 
the same way as the stellar disk: in an exponential disk 
with constant scale height equal to 6% of the disk scale 
length. A cuspy stellar bulge was included, with a mass 
1/3 that of the stellar disk. 

The mergers were distributed evenly across the 
phase space of initial geometries, with 2 simu- 
lations each from dircct+dircct, direct+retrograde, 
retrograde+direct, and retrograde+retrograde. (Di- 
rect and retrograde refer to the direction of rota- 
tion of each disk with respect to the encounter or- 
bit.) For each of these four combinations, we chose 
two fairly close encounters with initial pericenters around 
Rpericenter / Rdisk ~ 1 and 2. Each simulation is labeled 
self-evidently (e.g., ddl = direct— direct for disks 1 and 2, 
respectively, and R p /Rd ~ 1). Exact values of R per i, in- 
clination, and the direction of the line of nodes were cho- 
sen randomly within the constraints of direct/retrograde. 

With this set of conditions, we are not probing all of 
parameter space, but focusing on the case of a strong 
interaction between two present-day spirals of roughly 
equal mass. This is most relevant to the formation of 
ULIRGs and QSOs in the local universe (2J. Future 
work will focus on a wider range of parameter space. 

To get a sense of length scales, we note that for the 
best-fi t simulation to NGC 4676 considered in IBarnesI 
(2004), the linear scale was such that Rd = 3 kpc. This 
is comparable to the disk scale lengths o f spirals, which 
average several kpc ([Zaritskv et al.lll994l) . 

We calculate quantities of interest for each disk, as- 
suming that little mass transfer occurs from disk to disk. 
This is a valid assumption up to the second pericenter, 
when the disks begin to overlap appreciably. We thus 
constrain our analysis to before second pericenter. 

We include chemical evolution in our analysis by as- 
signing metallicities to the gas particles at the begin- 



ning of the simulation using the average gradient (pe r 
disk scale length) measured by IZaritskv et al.l (119941) . 
—0.2 dex/Rd- The dispersion in the IZaritskv et al.l 
(|1994D sample is 0.1 dex/Rd- We also add to each par- 
ticle a normally-distributed random component to the 
metallicity, with a = 0.1 dex and mean of 0. We then 
fix the metallicity of each gas particle for the duration of 
the simulation and follow the gas particles with time. 

We neglect enrichment from star formation, since we 
are focused on spiral galaxies in the nearby universe. 
Gas mas s fractions are sm all in local, massive spirals 
(~ 10%; iKannappanl 120 04) . so further enrichment will 
be small compared to that from previous generations 
of stars. (See [JS for further discussion.) We are be- 
ginning a study of simulati ons involving star formation 
( Springcl & Hernquist 2003) , which will be considered in 
future work. 

To study the radial metallicity distribution, at each 
time step we average azimuthally in radial bins of size 
0.1 Rd- We compute the gradient using a linear fit to the 
average metallicity in each bin, weighting each radius 
equally and fitting from to 5 Rd- (If, as mentioned 
above, Rd ~ 3 kpc for the progenitor disks, then the 
upper radius considered is ~15 kpc). This is comparable 
to the gradient measurement an observer makes, with 
H II regions distributed roughly evenly in radial space. 

3. RESULTS 

The finding that motivated this study was the ob- 
servation that interacting galaxies have suppressed nu- 
clear metallicities com pared to isolated ga l axies of simi- 
lar luminosity or ma ss dKewlev et alJl2006t lEllison et all 
20"08t iRupke et al l l2008t iMichel-Dansac et all 120081: 
Peeples et al.ll2009D . In Figure [TJ we show that this ef- 
fect is present in the simulations. The metallicity is un- 
changed prior to first pericenter, but it drops dramati- 
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Fig. 2. — Change in nuclear metallicity vs. number of gas par- 
ticles within the same radius. Each track represents the time evo- 
lution of a particular disk (disk 1 in blue, disk 2 in red) for each 
simulation, from lower left to upper right. The strong correlation 
is unambiguous, demonstrating that the undcrabundance is caused 
by gas inflow. 

cally after first passage. By second pericenter, the change 
in metallicity ranges from -0.1 dex to -0.3 dex, with an 
average of -0.2 dex. The effect is more dramatic if the ini- 
tial gradient is steeper (see below). The change in nuclear 
metallicity is correlated with the amount of gas inflow to 
the central region (Figure [5]), confirming the suggestion 
that it is due to the dilution of the central metallicity as 
lower-metallicity gas flows inward from the outskirts of 
the system. 

Turning to the metallicity profile of the entire disk, we 
show in Figure [3] the radial metallicity profile vs. time. 
Strong evolution from the initial distribution is evident, 
with the disks converging to a shallow or flat distribution 
over several disk scale lengths between first and second 
pericenter. The profile is not at all times characterized 
perfectly by a straight line, and tends toward a steeper 
negative slope in the innermost regions. 

Radial mixing also results in an increase in the disper- 
sion of the metallicity of the gas particles. From an initial 
dispersion of 0.1 dex, the dispersion within R/Rd < 0.5 
rises to a value of 0.15 to 0.25 dex at second pericenter. 
At larger radii, the change in dispersion is more dra- 
matic: it jumps to a typical value of 0.2— 0.25 right after 
first pericenter. This dispersion contributes to the jagged 
radial profiles at large radius and late times, which is 
also caused by small number statistics and discrete tidal 
structures. 

Accompanying the decrease in nuclear metallicity (Fig- 
urc[TJ) is an increase in metallicity beyond R ~ 3 Rd (Fig- 
ure [21) • This results from tidal tails that have spun out 
from the interiors of the progenitor disks, carrying small 
amounts of more metal-enriched gas from the galaxy cen- 
ters into the outskirts. This removal of metals is insignif- 
icant compared to the effect of radial inflow, in terms 
of the gas mass involved. Thus, the gas removal from 
the inner regions has no effect on the central metallic- 
ity. However, because of the relatively small gas mass 
at large radii, mixing of small amounts of more enriched 



gas easily raises the gas metallicity in these regions. 

To quantify the profile evolution, in Figure U we show 
the gradient as a function of time. We observe the follow- 
ing about the gradient: (1) it changes little before first 
passage; (2) it drops precipitously after first passage (by 
0.1 to 0.15 dex/Rdisk)', and (3) it experiences mild evo- 
lution after this first drop. The degree of change varies 
little with encounter geometry, but there is evidence that 
the ddl/2 models experience the most change, while the 
rrl/2 models experience the least. 

How does the change in gradient relate to the change in 
nuclear metallicity, which is more easily measured? Our 
models predict that, in the case of a major merger, the 
two quantities are loosely correlated (Figure [5]) . How- 
ever, the change in gradient proceeds more quickly than 
the change in nuclear metallicity, with the latter being 
most dramatic after the gradient has already flattened 
substantially. This is due to the large gas mass of the 
central region, which requires significant inflow for its 
metallicity to be strongly affected. 

In Figure [5j we demonstrate that the flattening of the 
gradient is not a strong function of the initial gradient, 
by comparing initial gradients of —0.2 and —0.4 dex/Rd, 
which correspond to the average and th e steepest gra- 
dients observed by Zaritskv et al.l (Qj)94). The range of 
gradients near second pericenter is to —0.1 dex and 
—0.05 to —0.2 dex, respectively. However, the change 
in nuclear metallicity can be much higher if the initial 
gradient is steeper. More precisely, the final change in 
nuclear metallicity (in dex) is roughly equal to the initial 
gradient (in dex/Rd). 

4. COMPARISON TO OBSERVATIONS 

Our results on nuclear underabundances com- 
pare favorably with observations. Previous work 
on a wide variety of interacting systems shows 
that underabundan c es are typic a lly . 1 — 0.3 dex 
(IKewlev et al.l [2001 IRupke et al.l [200l Ellison et all 
I2008t IMichel-Dansac et al.ll2008D . although higher un- 
derabundances exis t in some massive, blue systems 
(jPeeples et al.l [20091) . Our simulated offsets are consis- 
tent with most of these results, though an unknown 
fraction of the observed interacting systems represent 
weak or minor mergers. However, we c annot explain 
the s trongest deviations of 0.6 — 0.8 dex (jPeeples et al.l 
2009) unless the initial gradients a re steeper than those 
measured by iZaritskv et al.l (|1994D or the progenitors of 
these interacting systems were already underenriched due 
to lack of previous star formation (consistent with their 
blue colors). 

The present work is most applicable to mergers of 
equal-mass galaxies; the only systems that we are sure 
meet this criteria are ULIRGs. The metallicity offsets in 
these systems aver age 0.1 dex and range up to 0.3 dex 
(Rup ke et al.1120081 ) 1 . A gain, these results are in excellent 
agreement with the simulations, especially if we assume 
either that the ini tial gradients were so mewhat shallower 
than predicted bv lZaritskv et al.l (|1994| ) or that there has 
been modest re-enrichment due to star formation. 

1 T hese numbers are smaller than presented in IRupke et al.l 
(2008), and result from an updated analysis made possible by a re- 
calibrated mass- metallicity rela tion using the [N II] /[O II] metal- 
licity diagnostic (Kewlcy & Ellison 2008). This diagnostic is the 
most robust strong-line diagnostic (Kewlcy & Dopita 2002). 
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Fig. 3. — Radial metallicity profile at up to 10 time steps prior to second pericenter or the end of the simulation, whichever comes first. 
Blue (orange) dashed lines indicates early (late) times. The profile at first (second) pericenter is shown as a black (red) solid line. The 
profile for disk 2 is offset by -0.5 dex for clarity. Strong evolution toward a shallow or flat metallicity profile is evident between first and 
second pericenter. 



There is no published data on the metallicity gradi- 
ents of interacting systems. A study of weakly inter- 
acting spirals found no differenc e in metallicity grad i- 
ent compared to isolated spirals (jMarquez et al.1 [20031. 
but the [N II] /Ha metallicity diagno stics used in this 
study have significant uncertainties (jKewlev fc Dopital 
2002). We have be gun a program to observe the gradients 
in strongly interacting galaxies, to understand how the 
metal redistribution occurs over the course of a merger 
and to compare to simulations. 

It is clear that any significant galaxy-galaxy interac- 
tion will result in some change in the metallicity distri- 
bution of the galaxies involved. These effects may be 
important when interpreting the evolution of the mass- 
metallicity relation. At higher redshifts, when the merger 
and interaction rates w ere higher than at present (e.g., 
iKartaltepe et al.l [20071. the dilution of nuclear abun- 
dances may be a more ubiquitous phenomenon. The ef- 
fect will be especially significant if a population chosen to 
trace the evolution of the mass-metallicity relation con- 
sists of significant numbers of interacting systems. For 
instance, low-redshift Lyman Break Galaxy analogs show 
lower th an expected metallicity and evidence of a major 
merger (jHoopes et al.|[2007l lOverzier et al J 120081) . LBGs 
have been used to tra ce metallicity evolution at high red- 
shift (jErb et al.ll2006| ). However, higher redshift galaxies 
also have higher gas mass fractions, while our simula- 
tions assume a low gas mass fraction. Further work over 
a larger range of initial conditions is needed. 

One caveat to our results is the presence of ongoing star 
formation, which we have ignored. Early in the merger, 
this star formation occurs across the disk of the galaxy, 
and hence is unlikely to change our conclusions about 
gradients. Late in the merger, star formation is heavily 
concentrated toward the galactic center, and can raise 
the heavy element content in these regions. 

The magnitude of this effect depends on the initial 



metallicity and the fraction of the remaining gas that 
is consumed in star formation. Consider two spirals with 
initial n uclear metallicities of Z P) and a true oxygen yield 
of 0.005 (|Tremonti e t al. 200J, reduced by a factor of two 
to account for their high metallicity calibration). Since 
a significa nt amount of gas re mains at the late stages of 
a merger (lSanders et al.llf99l[ ). we assume that not more 
than half of the gas has been consumed prior to coa- 
lescence (more than this, and the initial gas mass frac- 
tions will h ave been unreasonably high) . In the closed 
box model (|Talbot fc ArnettJ [1971 . if 50% of the total 
gas mass is consumed in a solar-metallicity system, the 
metallicity will rise by at most 70%, or 0.2 dex. In this 
extreme scenario, enrichment by star formation could 
wipe out the nuclear metallicity underabundance pre- 
dicted by our simulations, and may affect the gradient 
as well, depending on where the star formation occurred. 
However, comparison to the data shows that ongoing en- 
richment is small up to the ULIRG phase, since under- 
abundances at the 0.1 — 0.3 dex level (or higher) are still 
observed in local interacting systems. This in turn im- 
plies that either we have assumed too much gas consump- 
tion in this example calculation (i.e., that the "extreme" 
case is not valid) , or that outflows (which are ubiquitous 
in these systems; iRupke et al.ll2005D have expelled metals 
produced in ongoing star formation, preventing further 
enrichment. 

Either way, our simulation results appear robust to on- 
going star formation. However, more sophisticated sim- 
ulations involving both star formation and metal redis- 
tribution by supernovae and outflows will be required to 
confirm this. We are beginning to analyze numerical sim- 
ulations of merging galaxies that include a prescription 
for chemi cal evolution due to star for mation and stellar 
feedback (|Springel fc Hernquistl 120031) . especially to un- 
derstand mergers at higher redshift with large gas mass 
fractions. 
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Fig. 4. — Radial metallicity gradient as a function of time. As suggested by Figure[3] the gradient abruptly turns shallow or flattens after 
first pericenter. The effect depends mildly on geometry, with retrograde+retrograde mergers showing the least effect. 



5. SUMMARY AND PROSPECTS 

We have presented numerical simulations that con- 
firm that nuclear metallicity underabundances observed 
in interacting disk galaxies are due to merger-driven 
inflow of low-metallicity gas from the outskirts of 
the progenitor disks. The magnitude of the effect 
is in agreement w it h what is observe d (~0.2 dex; 
Kewley etai] [2001 Hupke et al.l [200l Ellison et al.l 



20081 iMichel-Dansac et all 120081 : iPeeples et aiT T2009) . 
The nuclear metallicity drop traces closely the amount 
of merger-driven inflow, and is accompanied by an in- 
crease in the dispersion of gas cloud metallicities at a 
given radius. The metallicity drop and increase in dis- 
persion occur predominantly between first and second 
pericenter. 

Also evident between first and second pericenter is a 
dramatic flattening of the initial radial metallicity gra- 
dient. This flattening reflects the effects of gas redistri- 
bution over the galaxy disks. Such gradient changes will 
be observable in studies of H II region metallicities in 
strongly interacting systems. By comparing to detailed 
simulations, the distributions of gas-phase metallicities 
are certain to provide another way to constrain the evo- 



lutionary state of interacting systems. In particular, they 
are very sensitive to the gas motions that occur within 
the system. 

Our results apply primarily to local examples of major 
mergers of spiral galaxies. Analysis of simulations includ- 
ing more sophisticated chemical evolution prescriptions 
is currently underway, especially to probe a wider range 
of merger initial conditions (including mass ratio and gas 
mass fraction). It is important to understand the effect 
this could have on the determination of cosmic metallic- 
ity evolution, since galaxy mergers and interactions were 
more common at high redshift. Samples of high redshift 
galaxies that are dominated by interacting systems may 
bias metallicity determinations to low metallicities, com- 
pared to samples at low redshift dominated by isolated 
systems. Future work will explore this effect. 
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Fig. 5. — Change in nuclear metallicity vs. radial metallicity gradient. Each track represents the time evolution of a disk, from lower 
right to upper left. Blue (red) lines represent initial gradients of —0.2 (—0.4) dcx/iid- The two quantities are strongly linked, with some 
scatter. Steeper initial gradients result in lower nuclear metallicities, but still converge to shallow or flat gradients. 
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